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We show some evidences that the Southeastern Brazilian Continental Shelf (SBCS) has a devil's 
staircase structure, with a sequence of scarps and terraces with widths that obey fractal formation 
rules. Since the formation of these features are linked with the sea level variations, we say that the 
sea level changes in an organized pulsating way. Although the proposed approach was applied in 
a particular region of the Earth, it is suitable to be applied in an integrated way to other Shelves 
around the world, since the analyzes favor the revelation of the global sea level variations. 



I. INTRODUCTION 

During the late quaternary period, after the last glacial 
maximum (LGM), from 18 kyears ago till present, a 
global warming was responsible for the melting of the 
glaciers leading to a fast increase in the sea level. In ap- 
proximately 13 kyears, the sea level rised up to about 120 
meters, reaching the actual level. However, the sea level 
did not go up in a continuous fashion, but rather, it has 
evolved in a pulsatile way, leaving behind a signature of 
what actually happened, the Continental Shelf, i.e. the 
seafloor. 

Continental shelves are located at the boundary with 
the land so that they are shaped by both marine and ter- 
restrial processes. Sea-level oscillations incessantly trans- 
form terrestrial areas in marine environments and vice- 
versa, thus increasing the landscape complexity [l[. The 
presence of regions with abnormal slope as well as the 
presence of terraces on a Continental Shelf are indicators 
of sea level positions after the Last Glacial Maximum 
(LGM), when large ice sheets covered high latitudes of 
Europe and North America, and sea levels stood about 
120-130m lower than today [2j|. Geomorphic processes re- 
sponsible for the formation of these terraces and discon- 
tinuities on the bottom of the sea topography are linked 
to the coastal dynamics during eustatic processes asso- 
ciated with both erosional or depositional forcing (wave 
cut and wave built terraces respectively [H). 

The irregular distribution of such terraces and 
shoreface sediments is mainly controlled by the relation- 
ship between shelf paleo-physiography and changes on 
the sea level and sediment supply which reflect both 
global and local processes. Several works have dealt with 
mapping and modeling the distribution of shelf terraces 
in order to understand the environmental consequences 
of climate change and sea level variations after the LGM 

In this period of time the sea-level transgression was 
punctuated by at least six relatively short flooding events 
that collectively accounted for more than 90m of the 
120m rise. Most (but not all) of the floodings appear to 



correspond with paleoclimatic events recorded in Green- 
land and Antarctic ice-cores, indicative of the close cou- 
pling between rapid climate change, glacial melt, and cor- 
responding sea- level rise [3] • 

In this work, we analyze data from the Southeast- 
ern Brazilian Continental Shelf (SBCS) located in a 
typical sandy passive margin with the predominance 
of palimpsests sediments. The mean length is ap- 
proximately 250km and the shelfbreak is located at 
150m depth. It is a portion of a greater geomorpho- 
logic region of the southeastern Brazilian coast called 
Sao Paulo Bight, an arc-shaped part of the southeast- 
ern Brazilian margin. The geology and topography of 
the immersed area are very peculiar, represented by 
the Mesozoic/Cenozoic tectonic processes that generated 
the mountainous landscapes known as " Serra do Mar" . 
These landscapes (with mean altitudes of 800m) have a 
complex pattern that characterize the coastal morphol- 
ogy, and leads to several scarps intercalated with small 
coastal plains and pocket beaches. 

This particular characteristic determines the develop- 
ment of several small size fluvial basins and absence of 
major rivers conditioning low sediment input, what tends 
to preserve topographic signatures of the sea-level varia- 
tions. 

For the purpose of the present study, we select three 
parallel profiles acquired from echo-sounding surveys, 
since for all the considered profiles, the same similar se- 
ries of sequences of terraces were found. These profiles 
[1, S 03 are transversal to the coastline and the isobaths 
trend, and they extend from a 20m to a 120m depth. 

The importance of understanding the formation of 
these ridges is that it can tell us about the coastal mor- 
phodynamic conditions, inner shelf processes and about 
the characteristics of periods of the sea level regimes 
standstills (paleoshores). In particular, the widths of the 
terraces are related to the time the sea level " stabilized" . 
All this information is vital for the better understanding 
of the late quaternary climate changes dynamic. 

We find relations between the widths of the terraces 
that follow a self-affine pattern description. These rela- 
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tions are given by a mathematical model, which describes 
an order of appearance for the terraces. Our results sug- 
gest that this geomorphological structure for the terraces 
can be described by a devil's staircase [ll|, a staircase 
with infinitely many steps in between two steps. This 
property gives the name "devil" to the staircase, once an 
idealized being would take an infinite time to go from one 
step to another. So, the seafloor morphology is self-affinc 
(fractal structure) as reported in Ref. [H, , but ac- 
cording to our findings, it has a special kind of self-affinc 
structure, the devil's staircase structure. 

A devil's staircase as well as other self-affine structure 
are the response of an oscillatory system when excited 
by some external force. The presence of a step means 
that while varying some internal parameter, the system 
preserves some averaged regular behavior, a consequence 
of the stable frequency-locking regime between a natural 
frequency of the system and the frequency of the excita- 
tion. This staircase as well as other self-affine structures 
are characterized by the presence of steps whose widths 
are directly related to the rational ratio between the nat- 
ural frequency of the system and the frequency of the 
excitation. 

In a similar fashion, we associate the widths of the 
terraces with rational numbers that represent two hy- 
pothetical frequencies of oscillation which are assumed 
to exist in the system that creates the structure of the 
SBCS, here regarded as the sea level dynamics (SLD), 
also known as the sea level variations. Then, once these 
rational numbers are found, we show that the relative dis- 
tances between triples of terraces (associated with some 
hypothetical frequencies) follow similar scalings found in 
the relative distance between triples of plateaus (associ- 
ated with these same frequencies) observed in the devil's 
staircase. 

The seafloor true structure, apart from the dynamics 
that originated it, is also a very relevant issue, specially 
for practical applications. For example, one can measure 
the seafloor with one resolution and then reconstruct the 
rest based on some modeling |l4j. As we show in this 
work (Sec. |V|, a devil's staircase structure fits remark- 
ably well the experimental data. 

Our paper is organized as follows. In Sec. [Til we de- 
scribe the data to be analyzed. In Sec. IIII1 we describe 
which kind of dynamical systems can create a devil's 
staircase and how one can detect its presence in exper- 
imental data based on only a few observations. In Sec. 
IIV1 we show the evidences that led us to characterize the 
SBCS as a devil's staircase, and in Sec. |V]we show how to 
construct seafloor profiles based on the devil's staircase 
geometry. Finally, in Sec. IVI1 we present our conclusions, 
discussing also possible scenarios for the future of the sea 
level dynamics under the perspective of our findings. 
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FIG. 1: [Color online] Profiles (depth versus the distance to 
the cost) of the Southeastern Brazilian Continental Shelf. The 
arrows indicate the terraces considered in our analyzes. The 
profile shown with a thick black line is the profile chosen for 
our derivations, reproduced also in (B). The other two pro- 
files had their original position of the two axes shifted by a 
constant value such that one can also identify the terraces 
observed in the chosen profile in these other two. Note that 
a translation of the profiles by a constant value has no effect 
on any of the scalings observed. The reason of this mismatch 
between the profiles is due to the local geometry of the cost 
at the time the sea reached that level. 



II. DATA 

The data consists of the tree profiles given in Fig. [TJA- 
B). The profile considered for our analyzes is shown in 
Fig. [JJB), where we show the Continental Shelf of the 
State of Sao Paulo, in a transversal cut in the direction: 
inner shelf ("cost") — > shelfbreak ("open sea"). The hor- 
izontal axis represents the distance to the cost and the 
vertical axis, the sea level (depth), d. We are interested 
in the terraces widths and their respective depths. 

The profiles shown in Fig. [JJ were the result of a 
smoothing (filte ring ) process from the original data col- 
lected by Sonar [15| . The smoothing process is needed to 
eliminate from the measured data the influence of the os- 
cillations of the ship where the sonar is located and local 
oscillations on the sea floor probably due to the stream 
flows. 

Smaller topographic terraces could be smoothed or 
masked due to several processes such as: coastal dy- 
namic erosional during sea-level rising, Holocene sedi- 
ment cover, erosional processes associated with modern 
hidrodynamic pattern (geostrophic currents). For that 
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FIG. 2: Reconstruction of the terraces. The width of the ter- 
race is given by the minimal distance between the two parallel 
dashed lines. 



reason we only consider the largest ones, as the ones 
shown in Fig. [21 (located at d = —30.01m with the width 
of I = 6.06fcm). As one can see, the edges of the terraces 
are not so sharp as one would expect from a staircase 
plateau. Again, this is due to the action of the sea waves 
and stream flows throughout the time. To reconstruct 
what we believe to be the original terrace, we consider 
that its depth is given by the depth of the middle point, 
and its width is given by the minimal distance between 
two parallel lines placed along the scarps of the terrace 
edges. Using this procedure, we construct Table U with 
the largest and more relevant terraces found. 

We identify a certain terrace introducing a lower index 
ii in / and d, according to their chronological order of 
appearance. More recent appearance (closer to the cost, 
less deep) smaller is the index n. We consider the more 
recent data to have a zero distance from the cost, but 
in fact, this data is positioned at about 15km away from 
the shore, where the bottom of the sea is not affected by 
the turbulent zone caused by the break out of the waves. 
The profile of Fig. [2(B) was the one chosen among the 
other tree profiles because from it we could more clearly 
identify the largest number of relevant terraces [l5T | . 



n 


d n (m) 


l n (km) 


1 


-30.01 


6.06 ± 0.05 


2 


-41.86 


1.59 ± 0.05 


3 


-54.01 


2.93 ± 0.05 


4 


-61.14 


1.73 ± 0.05 


5 


-66.69 


2.21 ± 0.05 


6 


-74.33 


0.80 ± 0.1 


7 


-79.75 


0.80 ± 0.1 


8 


-85.30 


0.80 ± 0.1 



TABLE I: Terrace widths and depths. While the depths 
present no representative deviation, the deviation in the 
widths become larger for deeper terraces. The deviation in 
the widths is estimated by calculating the widths assuming 
many possible configurations between the placement of the 
two parallel lines used to calculate the widths. 



III. THE DEVIL'S STAIRCASE 

Frequency-locking is a resonant response occurring in 
systems of coupled oscillators or oscillators coupled to 
external forces. The first relevant report about this phe- 
nomenon was given by Christian Huygens in the 17'/i 
century. He observed that two clocks back to back in 
the wall, set initially with slightly different frequencies, 
would have their oscillations coupled by the energy trans- 
fer throughout the wall, and then they would eventually 
have their frequencies synchronized. Usually, we expect 
that an harmonic, P w i, of one oscillatory system locks 
with an harmonic, Q W 2, of the other oscillatory system, 
lcadingto a locked system working in the rational ratio 
P/Q m 

To understand what is the dynamics responsible for 
the onset of a frequency-locked oscillation, that is, the 
reasons for which a system either locks or unlocks, we 
present the simplest model one can come up with to de- 
scribe a more general oscillator. This model is described 
by an angle 6, which is changed (after one period) to the 
angle f(0). So, f{9) = 9 + Q. In order to introduce 
an external force in the oscillator modeling also possible 
physical interactions with other oscillators, a resonant 
term, g, is added into this model, resulting in the follow- 
ing model 



where 



f(0) = + n-g(0,K) (modi), 



g(e,K) = ^sm27r6. 



(1) 



(2) 



Despite this map simplicity, the same can not be said 
about its complexity [13]. Arnold (see ref. [HI) studied 
this map in detail aiming to understand how an oscil- 
latory system would undergo into periodic stable state 
when perturbed by an external perturbation. 

For K = 0, Eq. (p} represents a pure rotation, which 
is topological equivalent to a twice continuously differ- 
entiable, orienting preserving mapping of the circle onto 
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FIG. 3: A complete devil's staircase, obtained from Eq. |T}, 
for K=l. 



itself [Theorem of Denjoy, see Ref. [l9(]. Therefore, the 
simple Eq. ([T]) can be considered as a model for many 
types of oscillatory systems. In fact, Eq. ([1]) represents a 
more complicated system, a three-dimensional torus with 
frequencies Wi and u>2, when viewed by a Poincare map. 
Thus, f2 in Eq. ([1]) represents the ratio w±/w2- When 
W\/w2 = p/q (with p < q) is rational, this map has a 
period p motion and its trajectory, i.e. the value of 9, 
assumes the same value after q iterations. For K = 0, 
the so called winding number W is exactly equal to Q, 
i.e. W=p/q. 

For K ^ (non-linear case) W is defined by 



W(Q,K) = lim 

n — >oc 

where 



h(0 , K) + h(Gi,K) + ... + h(6 n -u K) 



h(9,K) = Q + g(6,K). 



(3) 



(4) 



For K < 1, Eq. ([T]) is monotonic and invertible. For 
K = 1, it develops a cubic inflection point at = 0. The 
map is still invertible but the inverse has a singularity. 
For K > 1 the map is non-invertible. 

Arnold wanted to understand how periodic oscillations 
would appear as one increases K from zero to positive 
values. He observed that a quasi-periodic oscillation, for 
an irrational Q and K = 0, would turn into a periodic 
oscillation as one varies K, from zero to positive values. 
He demonstrated that a periodic oscillation has proba- 
bility zero of being found for K = (rational numbers 
set is countable while the irrational numbers set is un- 
countable) and positive probability of being found for 
K > 0. He also observed that fixing a positive value K, 
the winding number W [Eq. ([3])] is a continuous but not 



differentiable function of fi, as one can see in Fig. [3l 
forming a stair-like structure. 

If W(f2, K) is rational, it can be represented by the 
ratio between two integer numbers as W — ^. At this 
point, the frequency £1 and the frequency of the function 
g are locked, producing the phenomenon of frequency- 
locking, when W(K, f2) does not change its value within 
an interval Afl (a plateau) of values of f2. In fact, smaller 
is the denominator Q, larger is the interval Af2. 

As one changes f2, plateaus for W rational appear fol- 
lowing a natural order described by the Farey mediant. 
Given two plateaus that represent a and a ^ wind- 
ing numbers, with plateau widths of Afli and re- 
spectively, there exists another plateau positioned at a 
winding number W within the interval [^-, ^] given by 



= P1+P3 
Q2 Q1 + Q3' 



(5) 



The Farey mediant gives the rational with the smallest 
integer denominator that is within the interval [^-, 

Therefore, the Ail(^) plateau is smaller than AQ(^J-) 
and Af2(-^|), but is bigger than any other possible 
plateau. Organizing the rationals according to the Farey 
mediant creates an hierarchical level of rationals, which 
are called Farey Tree. The plateaus Af2i and AO3 are re- 
garded as the parents, and ASI2 as the daughter plateau. 

The interesting case of Eq. ([T]) for our purpose, is 
exactly when K=l. For that case, one can find periodic 
orbits with any possible rational winding number, as one 
varies Q. What results in a zero (probability) measure of 
finding quasi-periodic oscillation in Eq. (H|) , by a random 
chosen of the fl value. Also, for K > 1, due to the overlap 
of resonances (periodic oscillation), chaos is possible. 

The devil's staircase can be fully characterized by 
the relations between the plateaus widths, and the re- 
lations between the gaps between two of them. While 
the plateau widths are linked to the probability one has 
to find periodic oscillations, the gaps widths between 
plateaus are linked to the probability one has to find 
quasi-periodic oscillations, in Eq. (JTJ) . 

There are many scaling laws relating the plateau 
widths [H, [2l| . There are local scalings, which re- 
late the widths of plateaus that appear close to a specific 
winding number, for example the famous golden mean 
Wg = 2~ 1 , However, we will focus our attention in the 
global scalings, which can be experimentally observed, 
and only for the case where K=l. For this case, we are in- 
terested in two scalings. The one that relates the plateau 
widths with the respective winding numbers in the form 
4j (the largest plateaus), and the one that describes the 
structure of the complementary set to the plateaus Aft, 
i.e. the structure of the gaps between plateaus. The 
structure of the plateaus is a Cantor set as well as the 
structure of the complementary set. Therefore, a charac- 
terization of these sets can be done in terms of the fractal 
dimension Dq 22] of the complementary set. 
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FIG. 4: Representation of the gaps S" , S' , and S, used to 
estimate the fractal dimension Do of the complementary set, 
using Eq. 0. 



The first scaling is [if: 



(7>3), 



(6) 



The second scaling relates the widths of the complemen- 
tary set as one goes to smaller and smaller scales. These 
widths are related to a power-scaling law whose coeffi- 
cient Do is the fractal dimension of the complementary 
set. For K = 1, we have that the fractal dimension of 
the complementary set is Do = 0.87. This is an univer- 
sal scaling. Since the complementary set of the plateaus 
represents the irrational rotations, the smaller is its frac- 
tal dimension, the smaller is the probability of finding 
quasi-periodic oscillation. 

For experimental data, the determination of Dq is dif- 
ficult to obtain because the dimension measures a mi- 
croscopic quantity of the plateaus widths, and in exper- 
imental data one can only observe the largest plateaus. 
Fortunately, an approximation D' for Dq can be obtained 
from the largest plateaus by using the idea proposed in 

M, 



the value of K through the largest denominator observed 

El 



Q> 



1 



(8) 



I -K 

Finally, we would like to stress that while in Eq. ([I]) the 
plateaus of the devil's staircase are positioned at wind- 
ing numbers defined by Eq. ([3]), in nature, devil's stair- 
cases have plateaus positioned at some accessible mea- 
surement. 

In the driven Rayleigh-Benard experiment [U, [25|, 
convection rolls appear in a small brick-shaped cell filled 
with mercury, for a critical temperature difference be- 
tween the upper and lower plates. As one perturbs the 
cell by a constant external magnetic field parallel to the 
axes of the rolls and by the introduction of an AC elec- 
trical current sheet pulsating with a frequency f e and 
amplitude B, a devil's staircase is found in the variable 
fi, the main frequency of the power spectra of the fluid 
velocity. As one varies the external frequency f e , stable 
oscillations take place at a frequency ratio f%/f e , for a 
given value of f e . In analogy to the devil's staircase of 
Eq. (P), f e should be thought as playing the same role 
of fl in Eq. |T]), and the ratio fi/f e as playing the same 
role of the winding number W. 

A devil's staircase can also be observed (see Ref. [26|]) 
in the amount of information H (topological entropy) 
that an unstable chaotic set has in terms of an interval 
of size e, used to create the set. To generate the unsta- 
ble chaotic set, we eliminate all possible trajectories of a 
stable chaotic set that visits this interval e. In analogy 
with the devil's staircase of the circle map, e should be 
related to Afi, while H to W. 

The first proof of a complete devil's staircase in a phys- 
ical model was given in Ref. 27], in the one-dimensional 
ising model with convex long-range antiferromagnetic in- 
teractions. 

In Ref. 28], it was found that a model for the El 
Nino, a phenomenon that is the result of a tropical ocean- 
atmosphere interaction when coupled nonlinearly with 
the Earth's annual cycle, could undergo a transition to 
chaos through a series of frequency-locked steps. The 
overlapping of these resonances, which are the steps of 
the devil's staircase, leads to the chaotic behavior. 



IV. A DEVIL'S STAIRCASE IN THE 
SOUTHERN BRAZILIAN CONTINENTAL SHELF 




where S', S, and S" are represented in Fig. [4] 

In case one has K = 1 (K < 1), we do not have a com- 
plete devil's staircase. In other words, winding numbers 
with denominators larger than a given Q are cut off from 
the Farey Tree. Using this information we can estimate 



The main premise that guides the application of the 
devil's staircase model to the Shelf is that the rules found 
in quantities related to the widths and depths of the ter- 
races obey the same rules found in a complete devil's 
staircase for the frequency-locked intervals Af2 and their 
rational winding number, W . Thus, we assume that the 
terrace widths l n play the same role of the frequency- 
locked intervals Af2, and the terrace depths play the same 
role of the rational winding number W . 



6 



In order to interpret the Shelf as a devil's staircase, 
we have to show that the terraces appear in positions 
which respect the Farey mediant, the rule that describes 
the winding number "positions" of the many plateaus. 
For that, we verify whether the position of the terraces 
at d n can be associated to hypothetical frequency ratios, 
denoted by w n = which respects the Farey mediant. 
In doing so, we want that the metric of three adjacent 
terraces respects the Farey mediant. In addition, we also 
assume that the larger terraces are the parents of the 
Farey Tree, while the smaller terraces between two larger 
ones are the daughters. Thus, for each triple of terrace, 
we want that 



d n +2 + d n W n+2 + W n 



d n -\ 



W n +1 



(9) 



One could have considered other ways to relate the 
depths and the frequency ratios. The one chosen in Eq. 
(J9j> is used in order to account for the fact that while d n 
is negative w n is not. 

From Eq. ([9]) it becomes clear that for the further anal- 
ysis the depth of a particular terrace does not play a so 
important role as the ratio between the depths of triple of 
terraces that contains this particular terrace. These ra- 
tios may eliminate the influence of the local morphology 
and the influence of the local sea level dynamics into the 
formation of the Shelf. Therefore, the proposed quantity 
might be suitable for an integrate analysis of the different 
Shelfs all over the world, specially the ones affected by 
local geomorphological characteristics. 

From the Farey mediant, we have a way to obtain the 
frequency ratios associated to each terrace, 



w n+ i = 



p n + p n+2 



(10) 



<?n + <7n+2 

Therefore, combining Eq. @ and Eq. (TIT))) , we obtain 



Pn + Pn+2 
In + Qn+2 



Pn + 2 , Pr, 
E 



(11) 



which results in 



(Pn +p n +2){E - l)q n q n +2Pn = PnVl+2 + &Pn+2, (12) 

where E is defined by 



d n +2 + d n 
d n +i 



= E. 



(13) 



We do not expect to have Eq. (112)) satisfied. We 
only require that the difference between the left and right 
hand sides of this equation, regarded as Se, is the low- 
est possible, among all possible values for p rn and q m 



n 


Pn 


q n 


1 


1 


8 


2 


2 


17 


3 


1 


9 


4 






5 






6 


1 


17 


7 


2 


35 


8 


1 


18 



TABLE II: Integers associated with the n considered terraces, 
with n = 1, .... 8. 



(with m — n,ri + 2), for a given E, with the restric- 
tion that the considered largest terraces are related to 
largest plateaus of Eq. {!]), and thus p m +2=Pm and 
qm+2=qm + 1, and Se <C 1. Doing so, we find the ra- 
tionals associated with the terraces, which are shown in 
table [III The minimal value of Se, denoted by min[<5e], 
is mm[Se(d 1 ,d 2 ,d 3 )}=0. 032002, with Pl / qi =l/8 for the 
terrace 1, and ^3/93=1/9, for the terrace 3. We also find 
that min[<5e(<i 6 ,d 7 ,d 8 )]=0.002344, with p 6 /q 6 =l/17, for 
the terrace 6, and ps/qs=l/18, for the terrace 8. These 
minimal values can be seen in Figs. EKA-B), where we 
show the values of 5e(di, d 2 , d^) [in (A)] and the values 
of Se(de, dT, ds) [in (B)], for different values of p and q. 
Using bigger values for p has the only effect to increase 
the value of Se. 

We have not identified rationals that can be associated 
with the terraces 4 and 5, which means that for p and q 
within p = [1,50] and q = [1,400], we find that Se > 1. 
We have assumed that they could be either a daughter 
or a parent. 

From now on, when convenient, we will drop the index 
n and represent each terrace by the associated frequency 
ratio. So, the terrace 1, for n=l, is represented as the 
terrace with w = 1/8. 

Table [TT] can be represented in the form of the Farey 
Tree as shown in Fig. [5J The branch of rationals in 
the Farey Tree in the form 1/q belongs to the most 
stable branch, which means that the observed terraces 
should have the largest widths. We believe that the other 
less important branches of the complete devil's staircase 
present in the data were smoothed out by the action of 
the waves and the flow streams throughout the time, and 
at the present time cannot be observed. 

Notice that as the time goes by, the frequency ratios 
are increasing their absolute value, which means that if 
this tendency is preserved in the future, we should expect 
to see larger terraces. 

In the following, we will try to recover in the experi- 
mental profile, the universal scaling laws of Eqs. © and 
©. Regarding Eq. ©, we find that / scales as l/g _3 <3 °, 
as shown in Fig. [71 which is the expected global uni- 
versal scaling for a complete devil's staircase. Regarding 
Eq. ([7]), and calculating S', S", and S using the triple 
of terraces with widths l(w = 1/8), l(w = 2/17), and 
l(w = 1/9), as represented in Fig. |H we find D'=0.89. 
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FIG. 5: [Color online] (A) Values of Se(d 1 ,d 2 ,d 3 ) and (B) 
Se(de, dj, d%), for different values of p and q. 8e is the differ- 
ence between the left and the right hand sides of Eq. (|12[) . 



FIG. 7: Scaling between the 1/g-terrace widths and the value 
of q. 
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FIG. 6: Farey Tree representing the frequency ratios associ- 
ated with the major terraces. 



Using the triple of terraces (n=3,n=4,n=5), we find that 
D'=Q.87. Both results are very close from the universal 
fractal dimension Do = 0.87, found for a complete devil's 
staircase. 



V. FITTING THE SBCS 

Motivated by our previous results, we fit the observed 
Shelf as a complete devil's staircase, using Eq. No- 
tice that the only requirement for Eq. (JTJ) to generate 
a complete devil's staircase is that the function g has a 
cubic inflection point at the critical parameter K = 1. 
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FIG. 8: Magnifications of the small box of Fig. [3] showing the 
plateaus of the devil's staircase of Eq. {1} that appear for the 
same frequency ratios associated with the triple of terraces 
it)=(l/8,2/17,l/9), in (A), and ^=(1/17,2/35,1/18), in (B). 



Whether Eq. ([T]) is indeed an optimal modeling for the 
Shelf is beyond the scope of the present study. We only 
chose this map because it is a well known system and 
it captures most of the relevant characteristic a dynam- 
ical systems needs to fulfill in order to create a devil's 
staircase. 

We model the SBCS as a complete devil's staircase, but 
we rescale the winding number W into the observed ter- 
race depth. So, we transform the complete devil's stair- 
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Depth - d{m) Distance to cost 



FIG. 9: The function F, which is a linear relation be- 
tween the frequency ratios associated with the terraces and 
their depths for the triple of terraces ui=(l/8, 2/17,1/9) and 
to=(l/17,2/35,l/18). 

case of Fig. [5] as good as possible into the profile of Fig. 
QjB), by rescaling the vertical axis of the staircase in Fig. 

m 

We do that by first obtaining the function F (see Fig. 
[5]) whose application into the terrace depth d(w) gives 
the frequency ratio w n = — associated with the terrace. 

For the triple of terraces w=(l/8, 2/17, 1/9), we obtain 

F(d[km\) = 0.14219 + 0.00057853d[fcm], (14) 

and for the triple of terraces w=(l/17, 2/35, 1/18), we ob- 
tain 

F{d[km\) = 0.080941 + 0.00029786d[/cm]. (15) 

Therefore, we assume that, locally, the frequency ratios 
are linearly related to the depth of the terraces. 

Then, we rescale the vertical axis of the staircases in 
Figs. UtA-B) and calculate an equivalent depth, d, for 
the winding number W by using 

d = F-\W), (16) 

We also allow tiny adjustments in the axes for a best 
fitting. The result is shown in Fig. [TDK A) f° r the triple 
of terraces w=(l/8,2/17,l/9) and in Fig. \MB) for the 
triple of terraces w=(l/17, 2/35, 1/18). We see that lo- 
cally, for a short time interval, we can have a good 
agreement of the terrace widths and positions, with the 
rescaled devil's staircase. However, globally, the fitting 
in (A) does not do well, as it is to be expected since the 



FIG. 10: [Color online] (A) Rescaling of Fig. [^A) in black, 
showing that the devil's staircase fits well the terraces with 
w=(l/8, 2/17,1/9) of the profile of Fig. \V^B), in gray. (B) 
Rescaling of Fig. |8jB) in black, showing that the devil's 
staircase fits well the terraces with ^=(1/17,2/35, 1/18) of 
the profile of Fig. DJB), hr gray- 



function F is only locally well defined and it changes de- 
pending on the depths of the terraces. Notice however 
that this short time interval is not so short since the time 
interval correspondent to a triple of terraces is of the or- 
der of a few hundred years. 

The assumption made that K = 1 is also supported 
from Eq. ([8]). Using this equation, we can obtain an 
estimation of the maximum value of K from a terrace 
with a frequency ratio that has the largest denominator. 
In our case, we observed w — 2/35. Using Q=35 in Eq. 
(J8|), we obtain K < 0.97. 

In Fig. fTOT A). we see a 1/7 plateau positioned in the 
zero sea level. That is the current level. Thus, the model 
predicts that nowadays we should have a large terrace, 
which might imply in an average stabilization of the sea 
level for a large period of time. However, this predic- 
tion might not correspond to reality if the sea dynamics 
responsible for the creation of the observed Continental 
Shelf suffered structurally modifications. 



VI. CONCLUSIONS 

We have shown some experimental evidences that the 
Southern Brazilian Continental Shelf (SBCS) has a struc- 
ture similar to the devil's staircase. That means that the 
terraces found in the bottom of the sea are not randomly 
distributed but they occur following a dynamical rule. 
This finding lead us to model the SBCS as a complete 
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devil's staircase, in which, between two real terraces, we 
suppose an infinite number of virtual (smaller) ones. We 
do not find these later ones, either because they have 
been washed out by the stream flow or simply due to the 
fact that the time period in which the sea level dynamics 
(SLD) stayed locked was not sufficient to create a ter- 
race. By our hypothesis, the SLD creates a terrace if it is 
a dynamics in which two relevant frequencies are locked 
in a rational ratio. 

This special phase-locked dynamics possesses a critical 
characteristic: large changes in some parameter responsi- 
ble for a relevant natural frequency of the SLD might not 
destroy the phase-locked regime, which might imply that 
the averaged sea level would remain still. On the other 
hand, small changes in the parameter associated with 
an external forcing of the SLD could be catastrophic, 
inducing a chaotic SLD, what would mean a turbulent 
averaged sea level rising/regression. 

In order to interpret the Shelf as a devil's staircase, we 
have shown that the terraces appear in an organized way 
according to the Farey mediant, the rule that describes 
the way plateaus appear in the devil's staircase. That 
allow us to "name" each terrace depth, d n , by a ratio- 
nal number, w n , regarded as the hypothetical frequency 
ratio. Arguably, these ratios represent the ratio between 
real frequencies that are present in the SLD. It is not 
the scope of the present work to verify this hypothesis, 
however, one way to check if the hypothetical frequency 
ratios are more than just a mathematical artifact would 



be to check if the SLD has, nowadays, two relevant fre- 
quencies in a ratio 1/7, as predicted. 

The newly proposed approach to characterize the 
SBCS rely mainly on the ratios between terraces widths 
and between terraces depths. While single terrace widths 
and depths are strongly influenced by local properties of 
the costal morphology and the local sea level variations, 
the ratios between terrace widths and depths should be a 
strong indication of the global sea level variations. There- 
fore, the newly proposed approach has a general charac- 
ter and it seems to be appropriated as a tool of analysis 
to other Continental Shelves around the world. 

Reminding that the local morphology of the studied 
area, the " Serra do Mar" does not have a strong impact 
in the formation of the Shelf and assuming that the lo- 
cal SLD is not directed involved in the formation of the 
large terraces considered in our analyses, thus, our results 
should reflect mainly the action of the global SLD. 

If the characteristics observed locally in the Sao Paulo 
Bight indeed reflect the effect of the global SLD, then 
the global SLD might be a critical system. Hopefully 
the environmental changes caused by the modern men 
have not yet made any significant change in a relevant 
parameter of this global system. 
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